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Abstract 

Coupled pair of PT-symmetric square wells is studied as a prototype of a quantum 
system characterized by two manifestly non-Hermitian commuting observables. Via 
the diagonalization of our Hamiltonian H ^ W and spin-like observable Vt ^ we 
demonstrate that there exists a domain of couplings where both the respective sets of 
eigenvalues n = 0, 1, . . . (energies) and a = ±1 ("spin projections") remain real. 
In such a "measurable" regime the model acquires a consistent probabilistic inter- 
pretation mediated by our selection of one of many available interaction-dependent 
scalar products. 
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1 Introduction 



One of the keys to the proposal of PT-symmetric Quantum Mechanics (PTSQM) by 
Bender and Boettcher [1] hed in the reahty of the spectrum of the imaginary one- 
dimensional oscillator well V(^cubic){x) ~ ix^- Although the rigorous confirmation of 
that fundamental as well as phenomenologically welcome property of the model has 
been delivered a few years later [2] , the proof remains rather abstract and complicated 
[3]. For this reason, a lot of parallel attention has been paid to the other, exactly 
solvable non-Hermitian potentials with real spectra [4]. People studied partially 
solvable (often called quasi-exact) analytic alternatives to V(^cubic) i^) ~ [5] as well 
as non-analytic square-well potentials of similar type [6] and their singular point- 
interaction limits [7]. 

Solvable choices proved particularly suitable for illustrative purposes. Their study 
clarified that PTSQM formahsm may be understood as a very natural extension of 
Quantum Mechanics, not asking for any new formulation of the "first principles". 
For a review we may recommend the recent dedicated Workshops' proceedings [8]. 

We intend to broaden the scope of the current PT-symmetric models beyond 
their popular ordinary differential equation (ODE) framework. We feel motivated by 
the observation that the majority of existing applications of the innovative PTSQM 
formalism concerns systems characterized by a single physical observable. We intend 
to fill the gap by an introduction of a model possessing a doublet of commuting inde- 
pendent obscrvables (sect. 2). Our square-well-type model is solvable and intuitively 
transparent (sect. 3). It exemplifies a number of generic features of PTSQM systems 
(cf. discussion in sect. 4). Its appeal and properties are summarized in sect. 5. 

A number of technical details is separated in Appendices A (an account of the 
perturbation representation of energies), B (summarizing the PTSQM formulism in 
a modified Dirac's notation), C (on norms) and D (on the single-channel projection). 



2 The model 

2.1 PT— symmetry and its generalizations 

The productivity of the counterintuitive PTSQM approach has been mainly revealed 
via studies of specific, concrete examples. Many of them proved too exceptional. 
Typically, one may recollect the elementary, exactly solvable spiked harmonic oscil- 
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lator of ref. [9] with its complete confluence of all the infinitely many 'exceptional 
points' defined as the couplings at which two neighboring real energies merge and 
complexify [10]. 

Among the less exceptional toy models, many conjectures have been deduced from 
discontinuous solvable potentials. The simplest, purely imaginary piece- wise constant 
potential of ref. [6] with single discontinuity contributed to our understanding of the 
mechanisms of stabilization of the real spectra [11]. Supcrsymmetric partners of this 
potential have been found obtainable by non-numerical means [12]. The study of its 
physical aspects and classical limit proved facilitated by its perturbative tractability 
[13]. A model-independence of most of these observations was confirmed by the 
long-range square- well model with two discontinuities [14], by the short-range model 
with three discontinuities [15] and by the harmonic oscillator decorated with two 
delta-function discontinuities [16]. 

All the above-mentioned solvable models offered an independent support for the 
inspiring conjecture of ref. [1] that the reality of spectra is related to the so called 
PT-symmetry of the potentials. This connects the observed absence of the complex 
energy eigenvalues with the invariance of the Hamiltonians with respect to the com- 
bined action of a complex conjugation T and a parity reversal P, VTH — HVT (cf. 
ref. [17]). 

The latter type of symmetry gave its name to all the PTSQM formahsm. Beyond 
the simplest ODE Hamiltonians, the emphasis on the parity and time reversal mean- 
ing of the operator VT may be weakened [18]. Thus, the parity V may be replaced 
by an arbitrary invertible self-adjoint operator P or, in a less confusing notation, 6. 
In parallel, the meaning of the complex conjugation T = T"^ may be extended to 
all the Hermitian-conjugation involutions AT'^ [19]. The PTSQM con- 

cepts become applicable to nonsymmetric operators and the PT-symmetry becomes 
re-interpreted as the property 

H^^eHe-\ e^e'' (i) 

called pseudo-Hermiticity of H [20]. 

An appreciation of the subtlety of the latter generalization requires a non-ODE 
model possessing more than one observable. We intend to describe here such a model 
in detail. 
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2.2 Two coupled channels 

The frequent use of the couphng of channels in physics [21, 22] attracted our attention 
to the partitioned 

with any parity-type invertible sub-operator Q = which is not necessarily involu- 
tive. 

A coupled pair of equal-mass particles moving in single spatial dimension inside 
a deep square- well box will be considered, exhibiting the symmetry (1) + (2) of their 
non-Hermitian Hamiltonian H = H(^kinetic) + H (^interaction)- In units h = 2m = 1 we 
shall have 

^ikineUO y ^ _^ j, ^irnteraction) ^ ^^^^^ ^^^^^ j • [^) 

The pseudo-metric (2) commutes with the kinetic (i.e., differential) operator HQ^inetic) 
so that the ^-pseudo- Hermiticity condition (1) will degenerate to an explicit defini- 
tion of H = Q~^V^Q and to the two ^— pseudo-Hermiticity relations 

wl = gWag-\ wl - gwj,g-\ 

Although we emphasized the generality of the operator g acting in the single-channel 
subspace, we shall simplify the discussion by a return to the common parity reversal 
in what follows, g — V. 

In order to select a specific channel-coupling interaction in (3) we shall pick up 
a maximally simplified, purely imaginary piece-wise-constant potential such that 
ReVa,b{x) = ReWa,b{x) = and 

Im Wa{x) = X > 0, Im Wbix) = Y >0, xE (-1, 0), 
lmWa{x)^-X, lmWb{x)^-Y, xe(0,l), 
ImVa{x) ^ImVb{x) ^ Z, x E (-1,0), 
Im Va{x) = Im Vb{x) = -Z, x e (0, 1) . 

This non-Hermitian model of the coupling of channels is defined in terms of its three 
real parameters X, Y and Z. This represents an immediate generalization of the 
single-channel square well of ref. [6] in which the spectrum happened to be real at 
all the not too large coupling constants [23]. Basically, we intend to prove the same 
for eq. (4). 
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3 Solutions 

3.1 Non-Hermit ian symmetry Q 

Hamiltonian H of eq. (3) enters the coupled-channel Schrodinger equation 

h( ^("^Ui^f ^("M. (5) 
V X{x) J \ X{x) J 

It specifies the bound states of the model when accompanied by the current asymp- 
totic boundary condition re-scaled to L = 1, 



X{x) 








(6) 



x=±L 

A merit of such a choice of the example is that its two-by-two Hamiltonian of eq. (3) 
commutes with the spin-like constant matrix 

^ / a;-M [X 
\uj J ^ Y 

It plays the role of another non-Hermitian "observable". Its eigenvalues a = ±1 are 
real and it exhibits also the 6'— pseudo-Hermiticity property, 

The existence of the symmetry Q implies that Schrodinger equation (5) may be 
complemented by the fixed-spin constraint 

n ( '^-^•"^ ]=a( '^-^•"^ ] . (7) 

V Xa{x) j \ Xc{.x) j 

It gives the relation between the channels at both the spin eigenvalues a — ±1, 

Xaix) = aujif^ix). (8) 
For the sake of brevity we shall mostly drop the subscripts ^ in what follows. 



3.2 Wavefunctions 

The connection (8) between the channels reduces the system of equations (5) into 
the single, a— dependent linear differential equation with the piece-wise constant 
coefficients, 

--r^(Pn{x) + [Va{x) + auWb{x)] (pn{x) = Enipn{x), X G (-1, 1). (9) 
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The necessary incorporation of the "asymptotic" boundary conditions (6) reduces 
further its general solutions to the ansatz 

1 C sin Kji{l — x), a; G (0, 1). 
The insertion of this ansatz in eq. (9) gives the linear relations 

E = kI + i{Z + a^/XY) = kI - i{Z + aVXY) (11) 

which define the energy and connect with at any given, fixed value of cr = ±1. 

Once we expect that the energies are observable we have to assume that all 
their values E — remain real. Vice versa, we are persuaded that for a certain 
fairly broad class of the coupled-channel non-Hermitian interactions and equations 
on a finite interval the general rigorous proof of the existence of a non-empty physical 
domain of parameters (where the energies remain real) may be based on the straight- 
forward extension of the proof delivered by Langer and Tretter in the single-channel 
case [24]. For our present purposes we shall feel satisfied by the less ambitious ap- 
proach paralleling simply the single-channel construction of ref. [6]. 

Within our physical domain V of the real X, Y, Z and E the inspection of 
eqs. (11) and (10) reveals that we may set kl — k*j^ — k — s — it. At a fixed spin a 
this re-defines 

E = s'^-f, Z = 2st-aVXY (12) 

in terms of some new pair of real parameters s and t. Wc may keep one of them 
positive (say, s > 0) while the second one lies on a branch of a hyperbolic curve, 

t^Uis) = y^Z,ff{a), Zeffia)^ Z + aVXY, a ^ ±1. (13) 

We see that the energies remain non-degenerate with respect to the spin cr in general. 

The quantization will be mediated by the requirement of the continuity of the 
wave functions (f{x) and xi^) and of their first derivatives at a; = 0. These conditions 
degenerate to the single pair of complex equations 

AsiuK — CsmK*, Akcosk — —Ck* cosK* . 

The first item fixes the normalization [A = Csmn* / sinre) while the elimination of 
C gives the complex constraint Re (k~^ tanre) = 0. It is equivalent to the single, 
(7— independent real secular equation 

s sin 2s + i sinh 2t = 0. (14) 
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Our construction of bound states is completed. They are determined by formulae 
(10) and (12) while their free parameters s and t must be fixed by the pair of eqs. (13) 
and (14). A few comments on the practical numerical and perturbative evaluation 
of the roots (s„, tn) may be found in Appendix A. 

4 Interpretation of the solutions 

One does not leave the Standard Textbook Quantum Mechanics (STQM) whenever 
feeling satisfied by the Hilbert space 'H{j,hysicai) where the scalar product [i.e., metric 
operator in its definition {a ^h) (physical) = (fl|0|^)] is kept trivial, Q(stqm) = I- In 
contrast, one is allowed and advised to admit a nontrivial metric within PTSQM 
framework, Q{ptsqm} 7^ I [25]. 

In the latter setting it is important to keep in mind that the non-Hermiticity of the 
operators of observables (i.e., of H and Q in our present illustrative example) might 
lead to some confusion in the standard Dirac's 'bra-ket' notation. For this reason, 
the slightly modified 'brabra-ket' notation of ref. [22] is advocated and summarized 
here in Appendix B. 

4.1 Physical metric 

The extended flexibility of PTSQM formalism is compensated by the necessity of 
an explicit construction of a consistent physical metric operator Q(ptsqm) 7^ It 
must be Hermitian (0 = 0''") and positive definite (0 > 0). The new freedom 
broadens the class of the observables H = C^(pysQM)' J — 1? 2, . . . , M which must be 
quasi-Hermitian, i.e., by definition, Hermitian in our new metric, 

H^f^eHQ-K (15) 

According to the review paper [26] the introduction of the observables of this type 
may be made in a mathematically consistent as well as phenomenologically appealing 
manner. Prom a pragmatic point of view, it proved particularly productive in nuclear 
physics where M > 1 as a rule. 

Requirement (15) looks difficult to satisfy, especially in less elementary quantum 
systems. The key to the technical feasibility of the transition STQM — > PTSQM 
has been found in the existence of an indeterminate, auxiliary pseudo-metric P (re- 
named as ^, in our coupled-channel eq. (1), in an attempt to avoid its easy confusion 
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with a very similar symbol V for parity). This means that one should speak, strictly 
speaking, about a ^T— symmetric Quantum Mechanics throughout our text. 

One usually proceeds in an opposite direction, from a (preferably, very sim- 
ple) pseudo-metric to metric. Thus, a complete set of eigenstates of a given set 
of ^— pseudo-Hermitian observables is constructed in the first step, and the proof of 
the reality of the energies is then added as quite a difficult task. In this sense, our 
present solvable model may serve as a source of an insight in the properties of the 
wavefunctions (cf. sect. 3.2) as well as of the domain of the reality of the spectrum 
(cf. Appendix A). 

On this background let us now address the correct physical interpretation of 
the theory. It is to be achieved via a specification of the physical metric, to be 
constructed as a Hermitian and positive definite solution = 0''' > of eq. (15). Its 
knowledge will enable us to treat any quasi- Hermitian operator A with the property 

= QAQ^^ as an observable. 

Once we start from the spectral representation (32) of H (cf. Appendix B), we 
may recall the property H^Q — QH and infer that 

e>= E \F,t)) Sp,r,E,a{{E,a\. (16) 
E,a,F,T 

The choice of the expansion coefficients S must remain compatible with eq. (15) and 
with the symmetry Jl. This gives the conditions 

SE,a,F,r {E* - F) = 0, SE,c,F,r (<7 - t) - 0. 

The spectrum of energies is assumed real so that the off-diagonal part of the array 
S must vanish. We may replace the quadruple sum (16) by the double-sum ansatz 

e = ^ \E^,a)) Sr,A{En,'j\- (17) 

It represents the formal solution of eq. (15) and contains the infinite sequence of 
arbitrary coefficients Sn,a- Whenever they do not vanish, 5'^o- 7^ 0, the operator (17) 
is formally invertible. 

The necessary [26] Hermiticity of is guaranteed when all the parameters Sn, a remain 
real. The necessary positivity of (which means its tractability as a genuine physical 
metric) will be achieved whenever all the coefficients remain positive, 5"^^ > 0. 
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We see that the acceptable metric (17) compatible with all the standard physical 
requirements exists (at least in the formal sense) and is non-unique. One is allowed 
to impose some other mathematical or physical requirements [18, 27]. 

4.2 Normalization conventions and the norm 

In our bound-state solutions \En,(T) of eq. (10) we are free to use any complex 
'normalization' constants C = Cn,a- The same freedom of choice applies to another 
series of the 'normalization' constants which would appear in the similar formulae 
for the 'left' eigenkets ((£■„, a\. This is unaffected by the observation of Appendix B 
that the respective definitions (29) and (30) are connected by the pseudo-Hermiticity 
property (1). One can only conclude that in all the non-degenerate cases the following 
proportionality rule remains valid, 

\E,a) - 0-^ \E*,a*)). 

For E = E* and a = a* at Z < Zcru "we may treat the latter rule as a definition of 
the eigenkets |-, •)) e 7i up to a normalization, 

1^, a)) = 9 \E, a) ^(^^f""^) , E^E^,Ei,... = real, a = ±1. (18) 

The exphcit solution of the left eigenproblem (30) is made redundant but the freedom 
in the choice of a convenient relative normalization (RN) factor survives. 

The same factor emerges in the following formula for the overlaps between eigen- 
states, 

{{En', cr'\En, a) = S^^>Snn' {En, cr\ e\En,a) , n, n' = 0, 1, . . . , a, a' = ±1. 

(19) 

Before its deeper analysis one may consult Appendix C which shows how our coupled- 
channel Hilbert space Ti. may be partitioned into two single-channel subspaces He- 
This partitioning is prescribed there in such a way that our (free) choice of an over- 
all normalization coefficient entering the right eigenstate \En,cr) is unambiguously 
inherited by its single-channel components |</7„) via eq. (33). In parallel, an inde- 
pendent choice of the coefficient in each left eigenstate ((£'„, (t| is transferred to its 
single-channel components by eq. (34). 

On this background we must inter-relate our full-space and subspace RN conven- 
tions. The factor ^^^^'"""'^ introduced in the full space H may differ from its subspace 
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partner \^°p^J;°"°-^^ of eq. (35). By construction, fortunately, both these quantities hap- 
pen to coincide [cf. eq. (37) and the rest of Appendix C for more details of the proof]. 
As a consequence, the partitioning enables us to replace the non- vanishing overlaps 
in (19) by the much simpler matrix elements (36). 

An inspection of the formulae (19) and (36) reveals that the absolute value of 
the self-overlap {{E, a\E,a) may be re-scaled to one by an appropriate choice of the 
'normalization' constants C = Cn,a in eq. (10). We also restrict the RN factors 
by the similar condition while their sign remains free, g^Z^"'^ = ±1- The key 
consequence lies in the fact that the sign of all the non-vanishing self-overlaps (36) 
is fully controlled by the sign of our optional RN factor. In an opposite direction, we 
have a freedom to prescribe such a specific set f?^^^""''' = ±1 which guarantees the 
positivity of all the 'special' self-overlaps. Thus, in the light of eq. (36) we simply 
postulate {{En, a\En, a) ^special) > 0, i.e., 

a g^Z:r'Hn\V\n) > . (20) 

After one evaluates the matrix element, this equation defines the dynamically deter- 
mined "physical" RN factors which make our basis biorthonormal, 

{{En', Cr'\En,Cr)(^special) — ^aa'^nn' ■ 

The (positive-definite) "physical" norm 

II \^n, cr)\\ (physical) 

sJ{{E,a\E,a) 

(special) 

of our bound states is obtained as a byproduct. This definition may be extended to 
all the elements of H via completeness relations (31) [18, 20, 28]. 



4.3 Quasi-parity Q 

In the literature people call the rc-scalcd RN coefficients a "charge" [18] or "quasi- 
parity" [29]. It is important to notice that in our present example their explicit 
determination is not difficult since the matrix elements {n\V\n) may be evaluated in 
closed form. Due to the purely trigonometric character of the wavefunctions (10) we 
have 

{n\V\n) = / sin K.*{x-\-l) sin k{—x-\-1) dx-\- sin K.{—x-\-l)sin k*{x-\-1) dx = 
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= — sin 2s cosh2f — ^ cos 2s sinh2i . 
2s 2t 

The n— dependence of this element is particularly transparent at the higher excita- 
tions with large s = 0{n) and small \t\ = 0{l/n). Our exact formula degenerates 
to its leading-order estimate 

^ {n\V\n) = -cos2s + f-) = (-1)" cosg„ + O f-) = {-!)" + O (-\ . 



AA* \nj \nj \n 

Thus, the "charge" or "quasi-parity" is specified by the closed formula 

at the higher excitations. 

Quasi-parities may now be interpreted as eigenvalues of a certain operator Q, 

Q\E^,a) = \E^,a)gir^K (21) 

We insert eq. (21) in (18) and deduce that 

{{En':Cr'\En,a)(^special) = {En':(^'\0 Q\En,a), n, n' = 0, 1, . . . , (7,C7' = ±1. (22) 

This identification defines an overlap of two different vectors using a specific scalar 
product in H. In the light of sect. 4.1 such a particular product corresponds to a 

particular metric operator, 

O (^particular) Q. 

Combining this relation with eq. (17) we get 

^^r"'^ = .r, — ■ (23) 

W-'-'n ) (7 I JZjfi , u j (special) 

Vice versa, the violation of the one-to-one correspondence (23) between the metric 
and norm would require an "artificial" introduction of an n— and a— dependence 
into our definition (21) of the operator Q. In spite of some formal merits of such 
a step [13] we are persuaded that the related "anisotropy" of both the operators Q 
and Q{particuiar) could hardly find a natural physical foundation. 

4.4 The crossings and degeneracies of levels 

Although the Hamiltonian H and "spin" Q (and wavefunctions) of our model depend 
on three parameters, its energy spectrum itself feels merely the influence of Z and 
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of the product XY . In the hght of eq. (28) the two sources of non-Hermiticity are 
the "internal" strength Z and the "couphng" strength \fXY . Still, the distinction 
between X and Y is nontrivial. At X — Y giving a "symmetrized" coupling of 
channels, the operator of symmetry Q becomes, incidentally, Hermitian. 

A strongly asymmetric decoupling of our model may be achieved by the two 
alternative hmiting transitions, viz., F 7^ 0, X ^ and X 7^ 0, y — > 0. Each of 
them suppresses just one of the channels [cf. (33)] . In both these limits the symmetry 
Q ceases to exist. One gets a;— >ooorl/a;— >oo and our present method of solution 
becomes inapplicable. 

A much more interesting limiting transition Z ^ (from both sides, i.e., Z ^ 
and Z 0^) converts our model into a coupled set of two Hermitian square wells. 
In this limit the violation of the Hermiticity of the whole system is merely caused 
by the channel-coupling terms. The energies E = — degenerate with respect 
to the spin since to-(s) = o"|to-(s)| so that, in fact, the neighboring a = ±1 levels 
cross at Z — 0. No point of the crossing is "exceptional" since the corresponding 
wavefunctions remain linearly independent. Their Wronskian W docs not vanish 
and both our observables H and fl remain diagonalizable. In contrast to some 
other solvable examples (say, to the harmonic oscillator of ref. [9]), no Jordan-block 
structures emerge in H. 

For the sufficiently small Zeff, all the similar observations may be made quanti- 
tative. Taking the ground-state n — and setting Zf,ff — 0{X) while Z — 0{X'^) for 
definiteness, we deduce that t — 0{X) [cf. eq. (13)]. Next we convert eq. (14) with 
s — 7r/2 + e and a small e — Eq [cf. eq. (24) in Appendix A] in the leading-order 

estimate of £ = At'^/n + All this transforms the definition of the energy into the 

following approximate formula 

As long as t^{s) = a \t(^{s)\, the ground state has the "spin" a = +1 at Z > and 
a — —1 at Z < while it becomes doubly degenerate at Z = 0. At this point the 
Wronskian easily evaluates to a nonvanishing constant, W = ACk* sm2K* . Hence, 
the two lowest states remain linearly independent at Z — 0. 
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5 Summary 



The appeal of virtually all the PTSQM constructions may be seen in a universality 
of their transition from a simple though indefinite pseudometric 9 to the correct 
physical and dynamically determined positive-definite metric ©. The procedure is 
counterintuitive and a number of open questions emerges. We designed our present 
less trivial coupled-channel example to clarify some mathematical subtleties (like 
the necessary conditions of the reality of the spectrum in non-Hermitian models), 
a deeper understanding of which requires, typically, a nontrivial application of the 
Krein-space theory [24] . 

We believe that the explanation of many interrelated subtleties of the PTSQM 
recipe may be facilitated via square-well models which offer one of the most econom- 
ical combinations of a transparent dynamical picture with an exact solvability of the 
underlying equations based on the usual matching technique. Our specific present 
example illustrates, first of all, a phenomenologically important situation where the 
dynamics is controlled by more observable quantities. 

An unexpected merit of our model has been found in a quick convergence of the 
auxiliary perturbation expansions of its cncrgy-lcvcl parameters s = Sn (such that 
E = — const/ s^) in the weak-coupling regime (i.e., for small strengths of the non- 
Hermiticity \Z\ and and/or in the quasi-classical regime (i.e., at the higher 

excitations with n ^ 1). Another, highly welcome byproduct of the square-well 
solvability emerged as a non-Hermitian spin-type symmetry Vt of H . It enabled us to 
reduce our nontrivial (viz., coupled-channel) Schrodingcr equation to its much more 
easily tractable "model-space" reduction. In parallel, the existence of the symmetry 
enabled us to analyze a level-degeneracy and level-crossing phenomena in a neat, 
non-numerical manner. 

The elementary algebraic structure of our model facilitated a clarification of one 
of the most puzzhng PTSQM requirements of keeping all the observables ©— quasi- 
Hermitian and, at the same time, ^— pseudo-Hermitian in the Hilbert space Ti.. The 
coupled-channel (i.e., partitioned) structure of the model enabled us to clarify the 
mechanism of this correspondence anew. In particular, we showed that the quasi- 
parity-based factorization © = as introduced in ref. [28] appears mathematically 
more natural than the alternative charge-based factorization © = CP of ref. [18], 
with P = ^ in our present notation. Indeed, while the quasi-parity is a symmetry of 
the Hamiltonian itself (we have [H, Q] = 0) , the formally equivalent introduction of 
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the charge C in [18] imphes that = 0. This means that the charge is merely 

a symmetry of an operator ^ H defined as a Hermitian-conjugate partner of the 
Hamiltonian. 

Due to the existence of the second, ^— pseudo-Hermitian and ©— quasi-Hermitian 
spin-hke observable Q in our model, another persuasive manifestation of a deep 
relevance of the symmetries of H has been revealed in the interrelations between the 
full space Hilbert space Ti. and its reduced, single-channel subspace Tic- Pars pro 
toto, the quasi-parity-related factorization Q = 9Q of the metric in full space TC has 
been proved accompanied by its analogue (40) using two relative quasi-parities 7^(ct) 
defined within the single-channel subspace TYc- 

The idea of the coupling of channels may turn attention to the systems treated 
perturbatively in more dimensions [30] as well as to non-pcrturbative explanations of 
the observed transitions between regular and chaotic classical and quantum motion 
controlled by the partial differential equations (PDE) [31]. Via our example, some 
existing confirmations of the internal consistency of the PTSQM theory may find 
their extension to the coupled-channel scenario. On this basis, "next" moves in the 
PTSQM development may be predicted as aiming at the non-separable PDE models 
[32] where some aspects of our model might inspire a more intensive exploration of 
the level-degeneracy patterns in non-Hermitian context [33] etc. 
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Appendices 



A. Perturbation series for the energies 

PTSQM models usually require some sufficiently efficient numerical description which 
is, typically, perturbative [34] . It may mediate an alternative or quicker insight even 
in the solvable models. Vice versa, the exact solvability of our present model offers 
an explicit verification of the approximative approaches. 

The predominantly trigonometric oscillatory character of the functions entering 
our secular eq. (14) enables us to locate and count all its physical roots, 

. = .„=(!i±i)£ + (-ir|, „ = 0,1,... (24) 

where the new parameter Qn remains small in the weak-coupling regime (i.e., for all 
the sufficiently small X, Y and Z) as well as at all the sufficiently large n. This 
enables us to abbreviate 

1 1 2Z,ff{a) 2Z,,ij(a) 

— Q=—, = a, rn = P = Q'g 



and to re- write our secular eq. (14) in terms of these new "small" parameters and a 
sign factor r = (—1)", 

Q = arcsin ( 2t sinh 2t] , 2t ^ . (25) 

As long as 

arcsin (x) — x -\- -x^ -\- -rr-a;^ -|- -rr-^x^ -\- . . . 

O "rrL/ X X ^ 

we may iterate eq. (25) and observe that the a— and /3— dependence of Q = <5(q;, /3) 
must acquire the following general asymptotic-series form 

oo 

g = g(a,/?) = a/?E(a,/?), S(a,/?)= ^ «''/3''cfe, (26) 

fc/=0 

where Cqo = 1- In practice, this series should be truncated in a as well as in (3 or, 
equivalently, in q. Nevertheless, as long as the size oi 13 = ga is dominated by a, it 
is sufficient to analyse this series as a power series expansion in the single "small" 
parameter a. It is also worth noting that at any fixed power of a there is always just 
a finite number of the related powers of q. 
For illustration, let us set 

E(a, 13) = 1 + cioa^ + Co 1 + c^o + Cn a^(3^ + Cq 2 + 0{a^) 
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and insert this formula in the re-arranged eq. (25), 

[l + T/32E(a,/3)]arcsinh|[l + r/32S(a,/3)]2i sin[a/3E(a,/3)]| = a. (27) 
As long as we can employ the regular Taylor series 

sm[aPE{a, P)] = a+(^cio + Cgi g^^ C20 + {cn - 1/6) + C02 a^-\-0 {a'^ 

the left-hand side of eq. (27) evaluates to a power series in our small parameters. 
The tedious though straightforward calculation converts the resulting equation into 
the infinite series dominated by the leading-order identity 

0= [-^ + + coi + ?>T q;^ + ... . 

It determines the first two coefficients, 

cio = ^, Coi = -3t. 
6 

Their insertion simplifies the next-order O (ct^) identity to the similar linear algebraic 
relation which defines the next set of the coefficients in S(q;, (3), 

In this manner one may continue the construction of the solution (26) to an arbitrary 
order in a. 

In the original notation we may now write down the second-order formula 

%Ztf. / 18(-l)"+i\ ^1// \ 



(n + 1)%3 3 + 1)67^5 \^ (n + 1)% V y(^ri + iyj 
etc. The convergence in l/(n-|- 1) is amazingly rapid and the role and weight of the 
non-Hermiticity decreases very quickly with the growth of the excitation n. 

In the single-channel limit where X — Y — Oit has been observed that the growth 
of |Z| = \Zeff\ makes some of the low- lying energies move towards each other. With 
the growth of the absolute value of Z they first pair (in fact, Eq and £"1) merges and 
complexifies beyond the critical value of Zcru ~ 4.48 [6, 23]. In the coupled-channel 
context we may repeat the same mathematical analysis leading, mutatis mutandis, 
to the conclusion that all the observable values of energies E = En and quasi-spins 
a remain real in the weakly non-Hermitian regime defined by the pair of inequalities 
\Zeff{a)\ < Zcrit , cr = ±1. They may be compressed into single condition 

I + \Z\< Z^rit ~ 4.48 . (28) 

In contrast to the single-channel case, the energy spectrum now ceases to be real at 
Z — ± (yZcrit — yfXY^i i-e., along two distinct surfaces in the space of parameters. 
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B. Modified Dirac's notation 



Our pair of operators H and samples a complete set of non-Hermitian commuting 
observables. These operators enter Schrodinger eq. (5) with quasi-spin constraint 
(7), i.e., in the Dirac's notation, the pair of equations 

H\E,a) = E\E,a), Q\E,a) ^ a\E,a) (29) 

As long as H and ^ O, the eigenkets |-, •) in (29) differirom the simultaneous 
eigenvectors of and In the spirit of ref. [22] let us now adapt the Dirac's 
notation to the non-Hermitian scenario and equip the latter elements of our Hilbert 
space 7i by the double delimiter. The apparently unmotivated complex conjugation 
of the energies and spins in their implicit definition 

\E, a)) = E* \E, a)), \E, a)) = a* \E, a)) 

becomes explained after the Hermitian conjugation which reveals the "left action" 
essence of these equations, 

{{E,a\H ^E{{E,a\, {{E,a\n ^ a {{E,a\ . (30) 

For this reason we shall prefer the use of the ket-vector form of the 'right' eigen- 
functions \E,a) in combination with the doubly delimited (or 'brabra- vector') form 
{{E, a\ of their 'left-eigenfunction' partners. 

Although both the latter sequences of vectors are defined, strictly speaking, in the 
two equivalent copies of the same Hilbert space of states H, our notation conventions 
will allow us to shorten the discussion here and there. As we already emphasized, our 
"redundant" version of the common Dirac's notation is transparent and proves more 
consistent in the non-Hermitian setting. Moreover, in the physical regime where 
E — E* and a — a* the two pairs of Schrodinger equations (29) and (30) imply the 
biorthogonality relations for their solutions which is easily written down now, 

{{E',a'\E,a){E' - E) =0, {{E',a'\E,a){a' - a) = 0. 

In the general non-degenerate case these rules only admit the non-vanishing over- 
laps at E' — E and a' — a. Vice versa, unless one of the self-overlaps vanishes 
accidentally, it is easy to derive the formal completeness relations 
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Their use enables us to treat our set of two sequences of states {{E^a\ and \E^a) 
as a biorthogonalized basis giving straightforward formal expansions of any element 
\a) = I ■ \a) en or \(5)) = l"^ ■ \(5)) en. Thus, one derives 

as two samples of an extension of the usual spectral representation to (arbitrary) 
operators emerging in the non-Hermitian coupled-channel context. These formulae 
will be needed in sect. 4.1. 



C. Partitioning of H into two subspaces He 

In eq. (19) we may employ the partitioned notation, 



K,a)=|''^"|' ^1, (7 = ±1, n = 0,l,... (33) 
\ifn) • aVX 

where the subkets |</?„) are cr— dependent solutions of eq. (9). In a left-action alter- 
native to this formula let us put 

((K,<t| = (<tVx((x„|, Vr((xn|), fT = ±l, n = 0,l,... (34) 

where the new subcomponents {{Xn\ are defined by a left-action version of eq. (9). 
More precisely, the left eigenstates {{Xn\ and the right eigenstates \(pn) correspond to 
the same reduced and spin-dependent parity-pseudo-Hermitian single- channel sub- 
Hamiltonian ^ 

H{a) = + K + (TioWt = V H\a) V 

which acts in the single-channel Hilbert subspace Tic and which is, in the language 
of ref . [1] , truly T'T— symmetric. 

The partitioning clarifies the structure of the spectral representations of the op- 
erators in our basis. All of them may be derived from the elementary projectors 



\E,a){{E,a\ 



\Vn) ■ X ■ {{Xn\ \Vn) " (^^fXY ■ {{Xr. 



entering, say, eq. (32). They may be understood as acting in two copies of Tic in 
n = nc®nc- We may abbreviate \ipn) = 1?^) and ((x„| = and collect the 
reduced Schrodinger equations, 

H{a)\n) = E,,\n), {{n\ H{a) = E^ {{n\, n = 0, 1, . . . . 



20 



They are to be solved in a single copy of TYc where their comparison leads to an 
alternative RN convention 

1^)) = V \n) aK°""'^ (35) 
paralleling cq. (18). We must check that and how both these normalizations remain 
mutually compatible. For this purpose we start from the RN definition (18) and add 
the partitioning (33) or, alternatively, start from the partitioning (34) and insert the 
definition (35) afterwards. In the former case we proceed via eq. (19) and get 

{{E,a\E,a) = 2a (^g^^'T ^ {n\V\n) . (36) 

In the latter case we have 

{{E,a\E,a) = 2a^/XY {{n\n) = 2a (A^^"^)* ^/XY {n\V\n) . 

A comparison of these two results reveals that 

\{cyptional) {optional) lo^\ 

^E„a = QE^a ■ \'^') 

Our two apparently independent RN constants must be chosen equal to each other. 
D. Quasi-parity in the subspaces Tic 

When we move to the single-channel subspacc Tic we encounter the two different 
bases { \n) } distinguished by the "external" parameter a. We have to fix a = +1 or 
(7 = — 1 in \n) = \ncr). This means that in a sub-space analogue of eq. (21) we have 
to define the two "reduced quasi-parities" as operators 7^(cr) in Tic with a manifest 
dependence on the spin, 

nia) \n) = nia) K) = \n) e^T""^ . 

At a fixed value of the spin a we obtain a subspace counterpart of eq. (22), 

((nV|na)(spe«aO = {n' „\Vn{a)\n^) , n, n' = 0, 1, . . . . (38) 

Although the spin- dependent product VTZ{a) plays just a not too important role of 
a subspace metric, its formal prolongation from Tic to the full space Ti is feasible and 
may be performed as follows. Firstly, one verifies that the action of the cr— dependent 
auxiliary operator 

S{a) = I I ^ann{a) 

aujlZ[a) 
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obeys the fixed-spin relation 

S{a)\Er,,a)^an{a)\n)n |^ ^ j = |K, a) ^^""'^ ■ 

A transition to the spin-independent formula will be then most naturally mediated 
by an introduction of the two two- by- two-matrix projectors Yl^j — {I -\- aQ)/2, 

s{a) n^lK, ^) = IK, ^) Qtr\ o = ±1. 

Wc may conclude that the spin-independent quasi-parity operator in the full space 7i 
may be defined by the formula Q = J2a=±i 5(cr) U^^. In the partitioned notation 
we may re-write this operator in the matrix form, 

^^1/ 7^(+l) + 7^(-l) o;-! [7^(+l)-7^(-l)] 
2 1^ u; [7^(+l) - 7^(-l)] 7^(+l) + 7^(-l) 

A return to another representation in terms of the projectors Ua is now possible, 

Q= Yl ^(^)na. (39) 

a=±l 

We see here that the two operators TZ{o') may be perceived as representing "reduced" 
quasi-parities in Tic- The new version of the factorization formula for the metric is 
delivered in the same spirit. 

In the light of eq. (39) this formula represents our factorized metric Q^spedai) = Qas 
a weighted sum of two factorized items equipped with the appropriate spin projectors. 



^{special) g 



iE (7 f ° , T)n,. (40, 

^ a=±i V 1 au ' J a=±i \ VTi{a) / 



This conclusion is compatible with formulae (36) and (22). 
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